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Abstract 

The dynamics of spin-lattice relaxation (sir) of metal-organic Pt(II) compounds 
is studied. Often, such systems are characterized by pronounced zero-field splittings 
(zfs) of the lowest-lying triplets. Previous expressions for the Orbach sir process do 
not allow to treat such splitting patterns properly. We discuss the behavior of a 
modified Orbach expression for a model system and present results of a fit of the 
temperature dependence of the spin-lattice relaxation rate of Pt(2-thpy)2 based on 
the modified expression. 
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Transition metal complexes with organic chelate ligands and their lowest ex- 
cited states are of potential use for solar energy conversion [1-7]. Recently, the 
processes of spin-lattice relaxation and the decay behavior of excited states 
have been studied experimentally for such systems in Shpol'skii matrices. [8- 
14] Of special importance are compounds with a Pt(II) central ion. Pt(II) 
systems exhibit many different types of low-lying excited triplets that in- 
clude metal-centered (MC) dd* states [15,16], metal-to-ligand-charge-transfer 
(MLCT) states [17-19], intra-ligand-charge-transfer (ILCT) states [10,13,20], 
ligand-ligand'-charge-transfer (LL'CT) states [6,7], and ligand-centered (LC) 
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Table 1 

Electronic origins E [ cm -1 ] (lowest triplet sublevel of Ti, lowest site), zero-field 
splittings [ cm -1 ] {AE ba : Energy difference between \b) and \a), AE^,: Energy dif- 
ference between |c) and \b)), spin-lattice relaxation times T g j r [ns] at 1.2 K, and 
transition types for various Pt(II) complexes with organic ligands 



Complex 


E 


AE ba 


AE cb 


T slr 


Type 


Ref. 


Pt(2-thpy) 2 a ) 


17156 


7 


9 


710 


LC/MLCT 


[9,21,14] 


Pt(2-thpy)(CO)Cl a ) 


18012 


0.055 


3.8 


3000 


LC/MLCT 


[14,24] 


Pt(phpy) 2 a ) 


19571 


6.9 


25.1 


390 


LC/MLCT 


[14] 


Pt(3-thpy) 2 a) 


18020 


13 


9 


» 25 


LC/MLCT 


[25,26] 


[Pt(bpy) 2 ](C10 4 )2 6) 


21237 


<1 


<1 


> 50 • 10 3 


LC/MC 


[23] 


Pt(qol) 2 a ) 


15426 


<1 


<1 


> 60 • 10 3 


ILCT 


[10,13] 


Pt(qtl) 2 «) 


13158 


<1 


<1 


>7000 


ILCT 


[13] 


Pt(phpy)(CO)Cl a ) 


20916 


<1 


6.4 




LC/MLCT 


[27] 


Pt(bhq) 2 c ) 


19814 


11 


28 




LC/MLCT 


[28] 


Pt(phpz) 2 a ) 


22952 


9 


7 




LC/MLCT 


[25] 


2-thpy : 2-(2-thienyl)pyridinate; 


phpy~: 


2,2'-phenylpyridinate; 3-thpy : 2- 



(3-thienyl)pyridinate; bpy: 2,2'-bipyridine; qol : 8-quinolinolate; qtl : 8- 
quinolienthiolate; bhq~: benzo[h]quinolinolate; phpz~: 2,2'-phenylpyrazinate. 

a ) In n-octane b ^ Neat material c ) In n-decane 

states with some MLCT and/or MC contribution [21-23]. In the following, we 
focus to Pt(II) systems with heterocyclic chelate ligands. 

As shown in Tab. 1, the low-lying triplets of these systems are characterized 
by a rather large variation of zero-field splittings (zfs) in the range from less 
than 0.1 cm -1 to about 40 cm -1 . The larger splittings are mainly due to 
spin-orbit coupling. For the same complex in different matrices, the lowest 
triplet states are shifted in energy (in many cases in the range of 200 - 400 
cm -1 ). The corresponding optical spectra show rich vibrational structure that 
may be well resolved (about 2cm _1 ) by choosing appropriate matrices and by 
employing methods of emission and/or excitation line narrowing. 

At low temperatures (several Kelvin), the processes of spin-lattice relaxation 
occurring between the triplet sublevels \a), \b), and |c) are relatively slow with 
relaxation times as long as hundreds of nano-seconds and even up to many 
micro-seconds (See Tab. 1 and Refs. [9-14]) due to the low density of phonon 
states corresponding to such zfs patterns. 

To discuss these processes, we assume that the perturbation V caused by the 
phonons couples the electronic states of the chromophore essentially linearly 
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Fig. 1. Processes of spin-lattice relaxation: a) Direct process, b) Orbach process, 
c) Raman process. 

(e.g. see Ref. [29, p. 228]) 

V = V^k + ... (1) 

k 



where is the strain corresponding to the phonon mode with wave vector k 
in the long wavelength limit. The matrix elements of Vy are denoted by V& Q = 
|(6|Vi|a)| and analogous expressions for V ca and V c b- The energy differences are 
AEba between \b) and \a), AE c b between |c) and \b), and AE ca between |c) and 
\a). The usual notation f3 = l/(ksT) for given temperature T and Boltzmann 
constant k B , and the abbreviations Cb a = C V b 2 a (AE ba f and analogous ones 
for C ca and C c b are also used. Here, the parameter (7 = 3/ (2nh 4 pv 5 ) is defined 
in terms of mass density p and (average) velocity v of sound of the matrix. 
The (AEba) 3 dependence of Cba should be kept in mind. 

The following relaxation processes (see Fig. 1) occur: 



Direct process: The rate is given by [30, p. 541], [29, p. 229] 



, (direct) _ , , 
h a,b — ^ab I ^ba 

= C ba coth((3AE ba /2) . (2) 

Here, k a b and k ba are the rate constants for the up and down processes, re- 
spectively, given by the expressions 



k a b — Cba 
kba = Cba 



1 

exp((3AE ba ) - 1 

exp(pAEbg) 
exp((3AE ba ) - 1 



(3) 
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Analogous expressions hold for the up and down rates fc& c , k c b, k ac , and k, 



Orbach process: The rate for this process vanishes for exponen- 
tially. It depends on the splitting pattern of the three involved states: If the 
energy separation AE ba of the two lower states \a) and \b) is much smaller 
than both the energy separations AE ca and AE c b to the upper state |c), then 
the well-known expression 

, (Orbach) _ ^C c bC ca 1 

a ' b " (C ca + C cb ) (e^-1) [V 



holds approximately for low T. This original Orbach expression is derived 
under the assumption that the energy differences are given by AE = AE ca = 
AE c b > 0. For a more general zfs pattern, the rate is given by the low-tempe- 
rature approximation [31] 

^(Orbach) k ac k c b ~\~ kbck ca kbckba 

k C a + k c b — kba 



with up and down rates as given in Eq. (3). The modified expression (5) 
contains Eq. (4) as a limiting case (see Ref. [31]). 

Raman process: For low temperature, the rate may be approximated by 

^(Raman) _ jjrpn ^ 



with a constant D and n = 5 for non-Kramers ions [32]. In the cases under 
study, this T 5 dependence fits the experimental observations [31] better than 
the T 7 dependence observed in other systems. 

The relative importance of the various sir processes is largely dependent on 
the size of the zfs and the energy separations to further electronic states. For 
instance, in systems like Pt(qol) 2 and Pt(qtl) 2 with a very small total zfs (see 
Tab. 1) and no further electronic states in the vicinity of Ti, direct and Orbach 
processes are expected to be very small due to the AE 3 dependence of these 
processes, and the Raman process is expected to dominate. Compare also Ref. 
[31]. 

The behavior of the above expressions is illustrated for a model system (with- 
out a Raman process) and with parameters AEb a = AE c b = 7 cm" 1 , v = 2000 
m/s, p = 1.1 g/cm 3 , Vbc = lOcur 1 , V ac = 20cm _1 , V a b = 3cm -1 . In Fig. 2, the 
relative errors of the approximations for both direct and Orbach process, i.e., 
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Fig. 2. Relative errors of the relaxation rate expressions with respect to Eq. (7) as 
a function of temperature T. Plotted are the errors of / c (* rec *) _j_ ^.(Orfeacft.) w -^.j 1 
(2) for the direct process in combination with the original Orbach expression (4) for 
different values of AE ( i) AE = AE cb = 7 cirT 1 , ii) AE = AE ca = 14 cm^ 1 , iii) 
AE = (AE ca + AE cb )/2 = 10.5 cm" 1 , iv) AE = AE &t = 5.4 cm^ 1 ) and v) with 
the modified expression (5). 

for the sum fc^ect) + k^ h rhacK) as obtained using Eq. (2) in combination either 
with Eq. (4) or Eq. (5), respectively. The errors are calculated with respect to 
the exact rate 



7 (Orbach+direct) 
K a,b 

= 2 i^ bc ^ ac ^ cb ^ ca ^ ba ^ a& ) 

~~ 2 i^^ bc ^ cb ~ ^ ab ~~ ^ ca ~~ ^ ac ^ ba ) 2 

\ 1/2 

+ ^(k c bk ca — k a bk ca — kb a k c b + kb a k a b) ) (7) 

for the three-level system that is obtained from the rate equations [31]. 

Applying the original Orbach expression, i.e., using Eq. (4) in combination 
with (2) for the direct process, the prefactor 2C c bC ca / {C ca + C c b) was com- 
puted from the model parameters, but different values of the parameter AE 
have been used: AE = AE c b corresponds to using the minimum distance of 
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state |c) to the states \a) and \b) (curve i) in Fig. 2), AE = AE ca corresponds 
to using the maximum distance (curve ii)), and AE = (AE ca + AE cb )/2 cor- 
responds to using the mean distance (curve hi)). The value AE = AE^ = 
5.4cm _1 is obtained by a least square fit of the exact data with one fit parame- 
ter AE (curve iv)), i.e., for the direct process and the prefactor of Eq. (4), the 
exact expressions were used during the fit. Interestingly, AE^ is less than any 
of the other differences of the energies. Alternatively, one could try to use the 
prefactor in Eq. (4) as an additional fit parameter. But then, one cannot hope 
to extract the model values of C cb and C ca from such a fit. Finally, curve v) in 
Fig. 2 was obtained using the modified expression (5) in combination with Eq. 
(2) for the direct process. Clearly, the modified approach yields much reduced 
errors over a large temperature range. Thus, Orbach's original expression (4) 
that was designed for a different pattern of the energy levels cannot be applied 
to a pattern with AE cb m AE ba for any reasonable choice of the parameter 
AE. 

We remark that similar results are also obtained for different choices of the 
parameters. For instance, for a value of v smaller by a factor /, the same 
results for the absolute rates would be obtained, if all the matrix elements 
of V 1 are also chosen smaller by a factor / 5//2 , e.g., for v = 1500 m/s and 
Vbc = 4.87 cm -1 , V ac = 9.74 cm" 1 , V ab = 1.46 cm -1 . Moreover, fixing all the 
other parameters, any rescaling of the three matrix elements by an arbitrary 
common positive factor yields the same error curves since we are dealing with 
relative errors and, under this scaling, all up and down rates k ab , k ba etc., and, 
hence, all sir rates in the model are multiplied by a common factor. 

It is of interest to present an example of the application of the above formalism 
to the spin-lattice relaxation observed for the lowest triplet of the Pt(2-thpy)2 
complex in an n-octane matrix. This compound is depicted in of Fig. 3, and 
some properties are collected in Table 1. The experimental spin-lattice relax- 
ation rate k^ slr ^ is obtained from the measured emission decay rate of state 
| b) by subtraction of the corresponding triplet deactivation rate to the ground 
state [31]. For the fit, we used Eq. (2) for the direct process, the modified 
expression (5) for the Orbach process, and Eq. (6) with n = 5 for the Raman 
process, i.e., for a T 5 low temperature dependence. As prefactor of the direct 
process, we used the low temperature limit of k^ slr \ The ratio of C ca /C cb can 
be obtained independently from time-resolved excitation spectra [9,31]. Also, 
all energy separations AE ba and AE cb are available from highly resolved spec- 
tra [9,21,14,31]. Thus, as fit parameter, only the prefactor D of the Raman 
process and the constant C ca remain. For such a two-parameter fit as displayed 
in Fig. 3, the result is highly satisfactory. 

A three-parameter fit based on the original Orbach expression (4) using the pa- 
rameters D, AE and the prefactor in Eq. (4) yields the value AE = 11.4cm _1 
(and a nearly doubled prefactor D for the Raman process in comparison to 
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Fig. 3. Fit of the spin-lattice relaxation rate k^ slr ^ as a function of temperature for 
Pt(2-thpy)2 in an n-octane matrix. Displayed are the contributions of the direct 
process (Eq. (2)), the Orbach process (using the modified expression (5)), and the 
Raman T 5 process (Eq. (6)). 

the fit displayed in Fig. 3). A similar value for AE was obtained in Ref. [9] by 
a somewhat different fitting procedure. Both these values are unphysical since 
they do not correspond to any of the observed energy differences (see Tab. 1). 
We remark that the present study was triggered by this difficulty of using the 
original Orbach expression (4). 

This result shows, as further ones presented in [33,31], that the use of the 
modified expression (5) for the Orbach process is necessary for a detailed un- 
derstanding of the dynamics of the spin-lattice relaxation for low-lying triplets 
of metal-organic transition metal compounds with their characteristic patterns 
of zero- field splitting. Thus, although the present study concentrated on Pt(II) 
compounds, the result should be applicable to a more general class of com- 
pounds, namely, the whole platinum metal group complexes (compare, e.g., 
the recent results [8,9] for [Ru(bpy) 3 ] 2+ ). 

Financial support by the Deutsche Forschungsgemeinschaft and the Fonds der 
Chemischen Industrie is gratefully acknowledged. 
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